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CL|' Abstract 

D , 

r^ • We propose and study at large N a new lattice gauge model , in which the Yang- 

Mills interaction is induced by the heavy scalar field in adjoint representation. At any 
dimension of space and any A^ the gauge fields can be integrated out yielding an effective 



\^ • field theory for the gauge invariant scalar field, corresponding to eigenvalues of the 



initial matrix field. This field develops the vacuum average, the fluctuations of which 
describe the elementary excitations of our gauge theory. At A^ = oo we find two phases 
of the model, with asymptotic freedom corresponding to the strong coupling phase (if 
there are no phase transitions at some critical N). We could not solve the model in 
this phase, but in the weak coupling phase we have derived exact nonlinear integral 
equations for the vacuum average and for the scalar excitation spectrum. Presumably 
the strong coupling equations can be derived by the same method. 



1 Introduction 

The complexity of QCD comes from the gauge field, which is neither completely random, 
nor completely classical. It is classical at small scales, according to asymptotic freedom, and 
random at large ones, according to quark confinement. When the large A^ approximation 
was invented in the seventies, we expected something to become classical at last, but such a 
master field was never found in spite of all efforts. At present the chance of solving large N 
QCD by a WKB approximation doesn't appear promising. 

Still, to our own surprise, we have found the master field. It becomes manifest, when we 
treat the gauge theory as induced quantum theory, in the spirit of Sakharov. The details 
are important here, it does matter what regularization to take, and what field to take for 
induction. 

We choose the scalar field in adjoint representation $(a;) whose eigenvalues (f)a{x) will 
later serve as a master field. As for the angular matrices Q{x) diagonalizing the matrix $(a;), 
those decouple after integrating out the gauge field. 

Our theory being nonperturbative, we have to specify the regularization. We have not 
invented anything new here, just taken the standard lattice theory, except we did not add 
the gauge field selfinteraction to our action 

S = J2 N tr iml^\x)- ^ $(a;)f/^(a;)<l>(a; + ^)f/^(a;) j ; tr $ = ; (1) 

X \ fi=l,2,..,D J 

The parameter rriQ is not quite the bare mass, but it could play the same role of adjusting 
the physical mass scale. In the next Section 2 we analyze the relation of this model to the 
conventional QCD, and obtain a scaling relation for the glueball mass scale as a function of 
mo. 

Then, in Section 3 we integrate out the gauge fields at every link, which yields effective 
interaction for the eigenvalues 4>a{x)- This interaction is essentially nonpolynomial, which is 
our only protection against the Gaussian fixed point for the master field. Investigating this 
interaction, we find a surprising relation between the well known Itzykson-Zuber integral 
over the unitary group and the two-matrix model, which describes the complete {p, q) table 
of conformal fields in 2D Quantum Gravity. Using modern large N technology, which is 
especially simple here, as there is no need for the double scaling limit, we find exact integral 
equation describing the weak coupling phase of the Itzykson-Zuber integral at given density 
of the eigenvalues of the scalar fields. 

The strong coupling phase is studied in Section 4. We could not find an exact analytic 
solution here, but rather studied the strong coupling expansion which we developed up to 
rriQ^^ . The general method for computation to any order of the strong coupling expansion 
is briefly described. 

In Section 5 we derive the classical equation for the vacuum average of the density of 
eigenvalues, and study it together with previous equations. As we find, one of equations can 
be solved exactly, which leaves us with only one nonlinear integral equation. We suggest 
that the critical phenomena arise when one of the endpoints of the spectrum goes to zero, so 
that two branch cuts of the weak coupling phase merge into one cut of the strong coupling 
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phase. 

In Section 6 we derive the effective Lagrangean for the scalar fields in the next large N 
approximation, corresponding to the Gaussian distribution for the vacuum fluctuations of 
eigenvalues. The corresponding quadratic functional diagonalizes in momentum space, but 
remains nonlocal in eigenvalue space. The mass spectrum is given by a certain linear integral 
wave equation, depending upon the master field. As one would expect here are in general 
infinitely many masses. Unfortunately though, in the weak coupling phase we studied, the 
hopping term in effective wave equation, producing the Laplace operator in continuum limit, 
vanishes in the large N limit by dynamical reasons. 

Of course, all these computations for the large N will have some physical meaning if there 
is no phase transitions happening on the way for some finite N, as it happens, in the Wilson 
type lattice QCD at A^ = 4. But even if this unwanted phase transition would appear in our 
model, there will be left still some possibilities in it to avoid it by an appropriate choice of 
the (nonlinear) potential for the field $. 

The remaining problems, and future directions are discussed in the last Section 7. We 
suggest numerical experiments to test its correspondence with QCD. 

2 Induced QCD 

Let us integrate out the scalar field Q in the functional integral of our lattice theory to induce 
the effective gauge field action 

JDUD^expi-S) <x JDUexpi-S^ndlU]) (2) 

The — expansion of the induced action 

mo '^ 

Sind[U] = -Tr In f S^^y - m^^ ^ Uf,{x) O Ul{x)6^+f,^y J (3) 

can be represented as a sum over lattice loops F of a scalar particle in external gauge field 
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where /[F] stands for the length of the loop, and f/[F] for the ordered loop product of U 
matrices. 

For weak smooth gauge fields the critical value of the bare mass would be at ml = d, but 
should we average over short wavelength fluctuations of the gauge fields, this critical value 
might shift to some other value ml. Near this value, we could take the continuum limit for 
the smooth part of the gauge field. 

Expanding in powers of the gauge potential, we would find the usual sum of one loop 
lattice Feynman graphs. It would be instructive to compute them explicitly, without first 

^In this Section we include the trace of $ as dynamic variable, which does not affect the leading order of 
the large N expansion. 



going to the local limit in the lattice action S and switching from the lattice to dimensional 
regularization. Let us outline this calculation. 

We represent f/^(a;) = exp(aA^(a;)) and choose the Schwinger gauge 

A^{x) = --F^^Xu + ... (5) 

where dots stand for the higher derivatives of the field strength. Note that the Abelian part 
of Afj_ drops from our effective action after taking absolute value of the loop product. So, 
this F^j^u is the SU{N) rather then U{N) field strength, as it should be. 
The leading terms in the adjoint loop product in the local limit 

|tr U[T]\^ ^N^- ^tr (F^,F„^) Y. ^^^Ay^ i^^ ' Vp) {^^ - Vu) (6) 

a;,j/cr 

where Ax = a;*^*-* — x*^*^^-* denotes the difference between two consequent points on the lattice 
loop. The identities Y^ Ax = were utilized in this relation. Now it is a straightforward 
exercise in lattice field theory to compute the sums over loops 



N 1 

-tr (F^^F^p) ^ ^^^^ ^;^ri Yl Ax^^Ay^ {xp - yp) {x^ - y^) (7) 

x,yCT 
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The sums over points x, y C F should be interchanged with the sums over lattice loops F, 
which yields the sum over all x, y on the lattice times the sum over numbers of links in the 
two parts TxyiTyx of the original loop F. As a result, the factor /[F] in the denominator 
cancels, and we arrive at independent sums over lattice paths Txy, '^yx- 

One could either go to momentum space and represent x — y as derivatives with respect 
to momenta, or stay in coordinate space, and use the asymptotic form of lattice propagator 

2«fr,.„l 2 r d'^P exp{tp{x-y)) r+°° d'^p exp{tp{x - y)) 



t. ° °^-(27r)^m2-E.>oCOs(p^) 7-oo {It^Y p' 



which is proportional to |x — yp '^. 

In the local limit we obtain the Yang-Mills action, with the bare coupling ^ given by 

the scalar one loop diagram with the ultraviolet cutoff A = a~^ , and the infrared one at the 
renormalized mass m^ = mg — m^. In four dimensions, with proper normalization 

1 AT , A2 

-^-— In— (9) 

gQ 9o7r^ m"^ 

The terms with higher powers of the Yang-Mills field strength would enter with inverse 
powers of m , but they would not depend upon the ultraviolet cutoff, as it follows from 
dimensional counting. This means that in the local limit A ^ oo one could safely neglect 
these higher terms. 

This calculation does not include the feedback of hard gluons. Now, once we obtained 
the induced coupling constant for the Yang-Mills theory, we could in principle estimate 



these effects. We should go baclcwards, from m to A, and take into account all the relevant 
interactions, including the 0^ induced by gluons. The resulting running gluon constant 
should blow up at A which would yield correct estimate for the relation between m and 
Qq. This relation depends upon unspecified bare 0*^ coupling, so that there seems to be no 
universality. 

We could have written from the very beginning the continuous version of our model which 
corresponds to the gauged scalar massive matrix field in the adjoint representation, without 
an explicit Yang-Mills term: 

L = ^tr ((a^^ + i [A^, $])' + m2$2 + Ao$') (10) 

By the naive change of the gauge we can gauge out the angular degrees of freedom of the 
$ field here. Than the gaussian integration over the gauge potentials leeds to an effective 
action for only the eigenvalues of $. But this action will suffer from the obvious drawbacks: 
different eigenvalues will be strongly attractive because of the Vandermonde determinants in 
the negative powers in the effective measure, which appear after the integration over vector 
potentials. Another problem is that from the very beginning we miss the diagonal degrees 
of freedom of the gauge field (it enters only in the commutator), which is of course bizarre. 
The resolution of all these problems will be found in the lattice version of the model where 
it will be obvious that some important contributions to the functional integral are missing 
in the continuous version, but are present on the lattice. 

Each of the terms in the lattice sum over loops has the same structure, as the usual Wilson 
action in adjoint representation. The fact that loops T are not elementary plaquettes, does 
not seem to be important, as this was the ambiguity of the Wilson theory anyway. 

It is important though, that these loops are small, which corresponds to heavy inducing 
fields. We are going to adjust the scalar mass so, that these loops would become much larger 
then the lattice spacing. In this case the effective action would be large, since the sum over 
loops would be close to divergency, which would provide us with nesessary large value of the 
bare lattice coupling /5 = ^ of induced gauge theory. On the other hand, the loops would 
be much smaller then the physical scale, so that their size would serve as an effective cutoff 
for induced QCD. 

One could go one step further in analogy with Wilson QCD. Namely, we could construct 
the Wilson-like strong coupling expansion for this multiloop gauge theory. The generic term 
of the Wilson strong coupling expansion corresponds to the closed surface made of elementary 
plaquettes. 

In our case these surfaces would be made from the variable loops, glued together side by 
side, in the same way, as elementary plaquettes of the Wilson expansion. Our surfaces are 
breathing, we have some "matter" on the world sheet of string. However, we do not want 
this matter to show up in the physical spectrum: we adjust the parameters so, that these 
scalar fields stay heavy. 

Coming back to the continuum limit, let us note, that the reason for leaving only the 
Yang-Mills term here repeats the original physical motivation for QCD to be the theory of 
hadronic interaction: regardless of what could be the quark and gluon constituents, the effec- 



tive theory above Grand Unification scales must be QCD as the only renormalizable theory, 
with the effective constant, running from something large at Planck scales to correct coupling 
at prehadronic scales. In our model, we take a particular choice of gluon constituents, which 
by universality and renormalizability must induce correct QCD. 

The scalar particle mass m would now serve as the ultraviolet cutoff for QCD, as it is 
the scale when the gluon interaction becomes nonlocal. Repeating the standard arguments 
of asymptotic freedom we would find the RG relation between the bare coupling and the 
glueball mass /i 

m^ 487r^ 

Comparing these two relations we find the scaling law: 

fi'^{ml-ml)^ (12) 

However, this relation was derived without taking into account the feedback of hard gluons, 
so that one could use it only as a hint. At the moment it seems, that there is no universal 
scaling index. The real issue is, of course, whether the model would escape the Gaussian 
fixed point at large scale. This would mean , that we are dealing with gauge theory. 

To summarize, this line of argument involved integration first over the scalar field to 
induce effective Yang-Mills theory, and then over gauge field to derive the RG law for the 
glueball spectrum. This is not the best way to solve our model, which was designed for the 
opposite integration order. 

3 Effective Interaction for Master Field 

Let us now integrate our theory in the opposite order, i.e. over the gauge fields first. This is 
possible in absence of direct gauge coupling in the bare Lagrangean. All we need , is the one 
link integral, which was computed first by Itzykson and Zuber[|l| (and was probably known 
before to mathematicians 0) 

/(^.x)=/z>f/exp(A'tr^f/xC/t)oc^^^|^^ep) (13) 

where 

A(0) = I[{^^ - ^j) (14) 

i<j 

is the Vandermonde determinant. 

Using eq. (/reflZ) we find from eq.(|lD the following partition function of our theory in 
terms of only eigenvalues (pi of the original scalar field: 

z^/nM~,wexp,-iv/2i:^w'^)A^ww)in^%g»^ (15) 

where < xy > denote the neighbouring vertices x and y on the lattice. 



Note thataboveinentioned naive consideration of the continuum model we would obtain 
only the diagonal term exp{N(j)i(f)j) from the determinant, with all the unwanted phenomena 
like the coUaps of the eigenvalues, whereas in the correct lattice version these problems are 
absent. 

Should we work with the physical A^ = 3 theory, this would be the explicit result for the 
effective interaction of the eigenvalues. However, in the large A^ limit, it is only half way to 
the result. 

The problem is obvious: there is the sign-alternating sum of A^! terms, each of the order 
of exp(A^^). They must cancel each other to some extent, to compensate the vanishing 
denominators at coinciding eigenvalues, since the integral of a positive finite function over a 
unitary group is always positive and finite. Instead of tracing these cancellations we use a 
different representation for the whole integral 



7(0, x) oc /■ d'^xd^'y ^ ^^^^''i'^uJ. .^) exp [n ^ x.y^ (16) 



A(x)A(y) 



where the integration contour over x{y) encircles the spectrum of eigenvalues (paiXa)- The 
proof is straightforward: one should compute the residues at the poles, which produce squares 
of Vandermonde determinants in denominator, combining with the same determinants in the 
numerator, after which the sign-alternating sum over all possible pole terms produces the 
determinant in the numerator. 

But this is the same as the two matrix integral 

y'rf^xd^yA(x)A(y)exp (NY.x,y, - Vl(x,) - V2{yi)\ (17) 

with singular potentials 

^iW = ^EM^-0a); v,{y) = ^j:Hy-Xa) (18) 

This integral can be computed at large A^ by means of the orthogonal polynomial 
technique 0. For the reader's convenience we reproduce this computation in Appendix. We 
assume here certain analytic properties of the potential, which we later find self consistent. 

The resulting equations of the two matrix models for the relevant case of x = 0? ^i = 
V2 = V read 

In /(0, (j)) = N\nhi + N'^ [ dt{l - t) In /(t) (19) 

Jo 

/^i = / 1^ / 1^ exp{Nxy - NV{x) - NV{y)) (20) 



dz V'{q{z ,t)) 

2m z^ 



m = -t^i^^^^^^^ (21) 



/ N 1 / f^^ '^"(g('2,^)) 



where the z integration goes counterclockwise surrounding all the singularities of the function 
V\q{z,t)). 

The last equation represents the integral equation for g(z, t) at fixed t, the second equation 
expresses t in terms of / , while the first one yields the answer for the large N limit of the 
Itzykson-Zuber integral. In the classical field equation of the next Section we shall need only 
the derivatives which we compute in the Appendix with the result 

I91n/(0,x)_ f'dt<f^—^— (23) 



N d(pa Jo J 2mz q{z,t) — 

In the usual case of the polynomial potentials the solution for q{z, t) has the form of a 
Laurent expansion 

q{z,t) = - + Y.q4t)z- (24) 

^ n 

which terminates at some integrals reduced to the residue at z = which produces a set of 
algebraic relations for the parameters qnit), f{t). 

The critical behaviour arises here by the standard mechanism, recently investigated in 
great detail in the matrix models of 2D Quantum Gravity [§], @], f^- Namely, the parameters 
must be adjusted so that the poles of the integrand in (^) pinch the integration contour 
at the symmetry points z = z~^. One could study, the requirement that the linear terms 
vanish in q{z, t) at z'^ = t = 1. 

^ = 0; .^ = 1. (25) 

In the vicinity of this point there would be powerlike singularities, the square root in the 
absence of further parameter adjustment. The structure of these singularities depends on 
the potential, which could be adjusted to produce any rational singularities. 

In our case the potential is not a polynomial. In general, we expect branch point singular- 
ities corresponding to finite support of eigenvalues. In this case the above integral equations 
should be investigated in the vicinity of these singularities together with the classical field 
equation for the density of these eigenvalues, which we derive in the next Section. 

Let us stress once again that the above formulas imply analyticity of the potential at 
the origin. One cannot apply them to the truncated expansions of the potential V'{q) in 
negative powers of q. The expansion in inverse powers of q corresponds to expansion in 
inverse powers of rriQ. This expansion is studied in the next Section. 

4 Strong Coupling Expansion 

We expect a phase transition in our model at a certain critical value rric of the bare mass 
mo, and we are interested in the strong coupling phase, when m,o \ rric. The above general 
solution applies to the opposite phase mo < rric, which will become clear below. In this 
Section we are going to study the strong coupling phase. Unfortunately, we could not find 
an explicit analytic solution in this phase, but we have found the strong coupling expansion. 



Let us rescale </) — i> — , y — ;> ^^ and expand our potentials in — 



mo ' '^ mo ^ ^ nil) 



V,(x) -. Inx - E ^^; \/2(:r) ^ Iny - J] ^^^ (26) 

where 

t„ = I rf0p(0)0'^; s„ = I rfxp(x)x" (27) 

are the moments of the distributions of eigenvalues. We shall assume that odd moments 
vanish, which corresponds to unbroken parity of our model. 

As it turns out, the above general solution of the two matrix model does not work in this 
expansion. In any order we find 

1 „/ N din 1(6, y) , , 

q{z) = -■ fit = 1 - t; -p^ = 28 

as exact solution. 

But we know that there should be — corrections. Straightforward expansion in powers 
of N tr (pUxU^^ in the original definition and U{N) group integration using Schwinger-Dyson 
equations yields, at infinite N, 

1 d In /^ ^ ^^^^ ^ ^^^3 ^ ^^^, ^ ^ , _,. ^29) 



N 



where 



Ts = -P+ ^ ^ '' ^ 30 

mo ^^0 



+ 



27 S2 ^2 - 30 S2 S4 t2^ + 7 S6 ta^ - 10 S2 ^4 + 10 S2 S4 t4 - 2 S6 ^4 



mo-*^^ 



^ _ -2S2^ + S4 -20 S2'^ t2 + 20 S2 S4 t2 - 4 S6 t2 

I4 — y \ 77 

mo' ?72o 



r. 



7S2 - 6S2S4 + S6 



rrio^^ 



The paradox is resolved by neglected singularities at small t in our formal solution. With 
a regular potential there are no singularities at small t so that the weak coupling phase 
is OK, but in the strong coupling phase the contributions from small t ~ -^ in the above 
integrals dominate. 

To see this, let us note that within the given order L of the — expansion one could write 

expW.,-V,W-Vi(,),) = ^^|H|||^^ (31) 

where 



are the L-degree polynomials in the inverse argument. The fact that there are higher powers 
of N in these polynomials does not bother us at the moment, since we are expanding in — 
at fixed A^, and we know that higher order terms would cancel in the final equations. 

This observation allows us to explicitely construct all the orthogonal polynomials with 
large order 

Piix) = x'-'R.ix); Qiiy) = y'~'R2{y)] I > L (33) 

Indeed, these are polynomials with proper normalization at infinity, and the integrals in 
orthogonality relation 

{Pn{x)Qm{y)) "^"^ / ^ / ^ exp(iVxi/)x"-^-iy™-^-i = d^^^K' K = ^^^ (34) 

correspond to f{t) = 1 — t at Nt > L. One can readily check that these polynomials Pi{x) 
are orthogonal to all powers x'^ ,k < I — 1. The y integral reduces to residue at zero, after 
which the x integral reduces to residue at infinity, which vanish. The same is true for Qi{y), 
with the obvious change of variables. 

Therefore, the derivative of our integral reduces to the sum of L terms 

d\nI{<P,x) _^/ Pn{x)Qniy) \ , . 

d<Pa tl\{^-<Pa)hnl ^ ^ 

— expansion of these terms can be explicitely computed using recurrense equations for the 
polynomials, working down to the zeroth and first order polynomials which are trivial. We 
have checked that the correct terms of — expansion are reproduced this way. 

mo L' L' J 

One may wonder why the lowest order integrals , such as hi, cannot be computed by the 
saddle point approximation at large A^. The answer is interesting. Within the — expansion 
there are multiple degenerate saddle points, so that one has to average over them, which 
leads to cancellations. One could also reproduce correct answers this way with some effort. 
The resume is that — expansion can be carried through to high order, which could be used 
for numerical computations. 

5 Classical Field Equation 

We are approaching the most surprising part of our theory, namely, the classical field equa- 
tion. The fact is, the integration over link variables U^{x) eliminates the angular variables 
of the scalar field, leaving us with only A^ eigenvalues at each cite of the lattice. Some 15 
years ago, without all the frustrating experience of large A^ QCD, one would not hesitate 
to apply the WKB approximation to the remaining functional integral over the eigenvalues. 
But now, the sophisticated reader might wonder where are the planar graphs. 

The answer is, we have already summed them up! The angular variables Q{x) of the 
scalar field $(a;) = Q^^{x)(f){x)Q{x) which we integrated at each site as the gauge part of 
the link integration, U^{x) -^ Q~^{x)U^{x)Q{x + /i) represent the majority of scalar degrees 
of freedom, responsible for the planar graphs. 

9 



For example, take the box diagram, corresponding to the one plaquette loop of the scalar 
particle. In conventional approach, one would expand in the hopping term at the original 
Action, and get the integral 



/Df/D^expf-^ATtr (m^$2(a:)) J []tr (<l>(i)f/(i)$(i + l)f/^(i)) ; $(5) = $(1) (36) 



i=l 



Integrating over $ first, one would get the usual Feynman graphs in an external gauge field, 
which then would be integrated out in the strong coupling approximation ( this is because 
we have not yet sumed up the scalar loops to find finite gauge coupling ). Integrating the 
box diagram in reverse order, we get a product of traces 



4 4 

i=l i=l 



(37) 



which yields zero since our scalar field is traceless. One could check that the same result 
comes about from the Feynman graphs of the scalar field. In a less trivial example of the 
double box, the integrations of 



4 



(3^ 



already yield a nonzero result. As in the previous example, the integration over link matrices 
yields the product of the traces in all vertices, but this time there would be traces of <l>^(i). 
At large A^ these traces can be computed in the WKB approximation for the eigenvalues (pa-, 
which for the Gaussian field corresponds to the famous semicircle distribution 



P(0) oc v^2 - mg02 



(39) 



Analyzing such examples, we come to the conclusion that link integrations take care 
of all the relevant quantum fluctuations at large A^, the scalar field as well as the gluon 
field. These are quite nontrivial integrals. They can be regarded as sums over A^! one loop 
WKB integrals around nontrivial "gauge vacua " , each integral producing two Vandermonde 
determinants in the denominator. All the "higher loop corrections " identically vanish, as 
a result of deep mathematical theorems about group manifolds. This means that we have 
summed over infinite number of "instantons" of our lattice theory, corresponding to one of 
each of A^! classical solutions at each link, plus the harmonic quantum fiuctuations. 

Let us now write down the classical equation for the eigenvalues of the scalar field, assum- 
ing that it is spatially uniform to fulfill translational invariance of the vacuum. Eliminating 
the angular variables, we get the square of the Vandermonde determinant at each vertex. 
The resulting saddle point equation for the constant master field reads 



2E' 



1 



2Nml^ 



-2d 



<91n/(0,x) 



(40) 



x=<t> 
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Let us consider the weak coupling phase first. Using formulas of the previous Section we 
find here at N = oo: 

i (v'ir) + nr)) = ™?^ - ij^ */^^ ; . - .(^.0 (41) 

where 0^ correspond to the opposite sides of the branch cut in the complex plane, corre- 
sponding to the spectrum of the eigenvalues. The equation holds only at this branch cut. 
Let us now introduce the analytic function -F(<^^) = (l)V'{(j)) (we assume that the parity 



transformation — >■ —0 of our theory is not spontaneously broken). The integral in (^T]) is, 
in fact, odd function of as g is an odd function of z for an odd potential, which allows us 
to rewrite the equation for F in manifestly symmetric form 

i(nn + nn)^mSf + <*-<*i'*/^^ (42) 

This is a boundary problem for the analytic function F{C,), with the conditions F{oo) = 1 
and F{0) = 0. The following Ansatz can be written down 

m = rnl{^-S{0) + d-dl'dtf^ f^f^^^^^^^^ ; S{^) ^ ^{^ - a-){r^ - b^) (43) 

This Ansatz satisfies the above equation, as S{C,^) = —S{^^) at the cut between a^ and 
6^. The fictitious singularities at ^ = g^ cancel between numerator and denominator of the 
integral. As for the conditions at and cxd, they yield the following equations 

^2 



m, 



0- .2 



'^"l/i^.w^) '"> 



a^ + 6^ \ Jo J 2mz S{q'^ 

The trivial solution, with a = or 6 = does not apply, since we assumed regular 
behaviour of F{C,) at the origin. Hence, we should choose finite solutions for a, b. This 
solution for the potential should now be substituted back into equations for q{z,t),f{t) of 
the previous Section. At present we do not know how to solve these equations, this seems 
to be a serious mathematical problem. 

As for the strong coupling phase, here we note that within the — expansion the master 
field equation coinsides with that of the planar limit of the one matrix model with f3 = Nd 
and effective potential Ueff{<p) such that 

f^e//(0) = (^ - ^2) - r40' - r60^ (46) 

Therefore, the — expansion is generated by the planar graphs of this model. The well 
known equations of the orthogonal polynomial the solution of spherical limit of the one 
matrix model can be applied now 

^ '-T2nRitr (47) 





11 



hn — S2n 



7d)2 Jo 



dt(mlR{t)Y (48) 



1 /•! e (4:R(t) - (f)"^) , , 

p(0) = - dt ) \[ ""J (49) 

27r ^0 J4i?(t) - 02 



The first relation provides an implicit equation for the function R{t), whereas the second 
equation yields the selfconsistency relations for the moments, involved in coefficients Tk 
defined in (0). 

These equations could be expanded in -^ 

„, , t (2n — 1)!! f X \ rriQ , -^ , , 

^ ^ = ;^^ + • • • ; S2n = Wttt + ■ • ■ ; p — = ir^ v/2^^ + . . . 50 

2mo (ra + 1)! Vmo/ 27r 

The eigenvalues are distributed from —h to 6 where }? = Amax R{t) ,0 < t < 1 within 
the — expansion. The phase transition to the above weak coupling phase would occur when 
the generating function 

F{i) = ^fdt ^ (51) 



li-AR{t) 

develops a singularity at ^ = 0. This singularity must come out from the second sheet, as 
there are no singularities at the first sheet within the — expansion. The density p(0) is a 
positive function, with a maximum at = 0, and therefore we expect p(0) to grow as some 
power of ttiq — ml. At the critical point the density diverges, after which we expect the gap 
from —a to a to arise in the eigenvalue distribution, in accordance with the weak coupling 
solution. 

In order to test this scenario one could expand the logarithmic derivative ^ in A- 
and Pade extrapolate to find the pole. The analytic computation of this series to high order 
is tractable, but exceeds our present scope. It seems to us that so far the basic physical 
properties of this remarkable model are more important then its numerical solution. 

6 Wave Equation 

Let us now discuss the excitations in our vacuum. In the harmonic approximation, corre- 
sponding to the next term of the large A^ expansion, we find in effective master field Action^ 

\ 1 ^ 

daab\ 5(t)a{x)5(l)b{x) - -riab5(t)a{x) Yl S4>b{x + fJ^) (52) 

/ fjL=—d 

where the matrix of second derivatives is computed in the Appendix 

= / ^ / ^ \ f53) 

/ 2Tn J 2711 {zi - Z2y{q{zi, 1) - 0a)(g(^2, 1) - 4>b) 

. . /(I 

Vab- 




f dzi r dz2 
J 2mf 2TnJT 



f{l)z^Z2y{q{z^,l) - <j)a){q{z2,l) - (Pb) 



^One may check that the Sab terms ah cancel. 
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As it is argued in the Appendix, the pole at zi = Z2 should be treated in the principal 
value sense, while the pole at f{l)ziZ2 = 1 lies outside the integration contours, which are 
little circles surrounding remaining singularities of the integrand. 

Note that the next corrections for the vacuum field do not display themselves in the 
leading order for the second variation matrix; it can be computed at A^ = cxd. The factor of 
A^ cancels when the field variation is rescaled Scj) = ^^^ and all the sums are converted to 
the integrals in the local limit 

^2 = fd'^xC (54) 

£ = I rf0p(0) J t/0'p(0')V^(0) {b{<P, 0')^(0') + C(0, 0')5V(0O 

where 

^ _ /(I) / / dzidz2 1 

2 ffi2myil-f{l)z^Z2nq{z,,l)-<P){q{z2,l)-<P') ^ ' 

B = , \^ + dC - d ^ ^ "^''"^'^ ' 



- 002 / / (27rz)2 (^1 - Z2)2(g(^i, 1) - </')(g(^2, 1) - </>') 

and double poles are to be treated in the principal value sence. 
The particle masses Mi are to be found from the wave equation 

d<P'p{<P') (C(0, cp')Mf - i?(0, 00) ^(00 = (56) 

which in general has a discrete spectrum, for the finite support of eigenvalues. Clearly, the 
bare mass mo should be adjusted to make the physical mass scale small in the lattice units 
we are using. Presumably, the scaling law of a type (1T2|) would come out of this equation, 
or from the similar equation of the strong coupling phase. The latter equation is yet to be 
derived, which is a challenge to large A^ experts. 

The serious problem with the weak coupling solution under consideration is that appar- 
ently /(I) = so that there is no hopping term in effective Action, and therefore no mass 
spectrum. 

The arguments for vanishing /(I) are very simple. As discussed in the Section 4, the 
orthogonal polynomials of high enough order can be explicitely constructed. In particular, 
for L > AT 

N N 

Pl{x) = x^-^-i n(^ - ^0; Qdv) = y''-''-' Uiv - y^) (57) 



i=l 1=1 



SO that 



]\tN-L , T • 



which provides us with boundary value /(I) = 0. 
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7 Discussion 

So, is this the long-anticipated exact solution of large A^ QCD? Not quite, but we are getting 
closer to it. It remains to be proven that this theory correctly induces QCD, which would 
require either an exact solution in the strong coupling phase, or numerical simulations. In 
our opinion, both are highly desirable. 

As for numerical simulations, these would be especially simple at iV = 2 where there are 
only two eigenvalues of opposite sign ±0, and the functional integral reduces to 

n /■=°#(x)0(x)=exp {^2ml^(xf) H ""'^tM'" f^^' 

X J^ ^ ^ <xy> m.x)<p{y) 

We do not see any reason why our model should fail to induce QCD at A^ = 2 as well as it 
does at A^ = oo. Clearly, the WKB approximation is not supposed to work so well here, but 
at least we could expect the vacuum average of 0. The restriction > could be taken into 
account by the change of variables = exp(A) where A varies from — cxd to +cxd. Still, the 
potential for the A field remains essentially nonlinear, so that it could avoid the zero charge 
problem. The model looks too simple, and maybe some flaw in it would be found, but it is 
worth a try. At least we would get some insight and see what could go wrong. 

Then one could try the SU{'i) model. What are the phases of this model? How does it 
compare with the usual QCD? Maybe the scaling laws would be easier to obtain numerically 
than the exponential laws of asymptotic freedom? The scaling indices in, say, the 3D Ising 
model are reproduced by state of the art lattice simulations with several digits, whereas 
asymptotic freedom has never been correctly reproduced, to the great embarrassment of the 
lattice gauge community. So, why not try our model instead? It does not look much harder 
than the Ising model. 

As for exact solution, the part of the theory which we miss badly is the strong coupling 
representation of the Itzykson-Zuber determinant in the large A^ limit. We feel that such 
a solution is just around the corner, but all efforts to find it from the two matrix model 
representation have failed. 

Still, we could approach the critical point from the weak coupling side, which we did 
in the previous Section. There are explicit nonlinear integral equations for the eigenvalue 
density and the wave equation for the glueball spectrum to solve, which is another interesting 
numerical problem. 
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A Itzykson-Zuber determinant via Orthogonal Poly- 
nomials in Two Matrix Model 

Let us reproduce and slightly generalize Mehta's solution of the two matrix model by means 
of orthogonal polynomials. The basic idea is to introduce orthogonal polynomials Pn{x), 
Qm{y) with properties 

{Pn{x)QUy)) = KSnn.; {F{x,y)) = j ^J ^^F {x , y) e^N xy - NV^{x) - NV^^y)) 

(60) 
The integration contour surrounds the singularities of the measure, which in our case are 
simple poles at x = 0^ and y = Xb- The normalization of polynomials is such that they tend 
to Pn{x) — > x"'~^ ,Qm{y) -^ y^-^ at infinity. Note that we shifted the index by one unit 
with respect to traditional definition - this simplifies formulas below. The Vandermonde 
determinants can be expressed as determinants of these polynomials 



A(a;) = deiP,{xi) = e{„} J]^™.!^*) (61) 

A(t/) = det Qj{xi) = e{m} TT QmXyi) (62) 

where e{„} denotes the unit antisymmetric tensor with A^ indices ni, . . .71^. 
With these substitutions our original integral reduces to 

TV 

I OC €{„}£{„} Y[ {PnXXi)QmXyi)) = N\\{hn (63) 

i n=l 

A convenient set of parameters, with smooth limit at large A^ is given by the coefficients 
qi(t) ,pi(t) of the expansion ( t = -^ ) 

xPn{x) = Pn+l{x) + Y.Pl(^)Pn-l{x)- yQn{y) = Qn+M + T.1l(t)Qn-liy) (64) 

I I 

which for finite A^ terminates aX I = N but in the large A^ limit becomes an infinite expansion. 
The equations for these parameters follow from identities 

= (P^,(x)g„(t/)) (65) 
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Integrating by parts in x we find 

{Pn-i{x){V;{x)~y)Qn{y)) = (66) 



Now, using the expansion (|J) together with orthogonahty condition (pOD we arrive at the 
relation 

q,it) = h-'_^{F,_iix)V;{x)QM) (67) 

and the similar relation for pi{t) 

Pi{t) = h-\ {Qn-i{y)Vi{y)P,,{x)) (68) 

There is one more relation, for the /i„ coefficients. It follows from the identity 

= ((P^+i(x)-nP„(x))Q„(y)) (69) 

after integration by parts of the first term 

= ((P„+i(x)iV(V;'(x) -y)- nPnix)) QM) (70) 

and using orthogonality relation 

K+^ = --K + {Pn+i{x)Vi{x)QM) (71) 

At this point it is convenient to introduce the ratio 

fit) = %^; ^ - ^ (72) 



'"n 



and the normalized P , Q polynomials 

Qn{y) = %^; Pnix) = ^ (73) 

{Pn{x)Qrn{y)) = {Q rn{y) Pn{x)) = 5nm (74) 

and rewrite above equations as follows 



N-l 



k 



lnJ = lnA^! + A^ln/ii+ ^(A^-A;)ln/ — (75) 

fc=i \^/ 

fit) = -t + {Pn+i{x)Vi{x)Qn{y)) (76) 

q,(t) = {Pn-i{x)Vi{x)QM) n / (^ - ^) (77) 

Pi{t) = {Qn-i{yK{y)Pn{x)) n / (t - ^) (78) 
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As for the derivative of the free energy, it reduces to the following 

d In /(0, x) _^ / Pn{x)Qn{y) \ , . 

Now, assuming existence of the large A^ limit of / we arrive at the asymptotic formula 
for our integral 

^^^ '^ dt{l-t)\nf{t) (80) 



'1) 



m JO 

The products in above equations converge to 



ri/(-a-/'<"(-^'"^)-^^-/'<" 



As for the normalized matrix elements, those could be computed at large A^ in the WKB 
approximation, using the shift operator 

. = exp (-1) (82) 

and neglecting its commutators with t 

[z,t] = -1^0 (83) 

The recurrent equations for our polynomials can be rewritten as operator relations 

yQn^kiy) = qiz, t)z^Qn{y)\ xPn-kix) = p{z, t)z^Pn{x) (84) 

with 

q{z,t) = - + Y.^i{t)z'- piz,t) = - + J2Mt>' (85) 

^ I ^ I 

Furthermore, 

V;ix)Pr,-lix) = J2 Cl^Pn-m{x) (86) 

m 

where the coefficients Cim can be computed by WKB formula 

Cim = f^/-'^-'VMz,t)) (87) 

and likewise for Q. The z integration contour surrounds all the singularities of the integrand, 
which for the polynomial potential would reduce to the residue at z = or, equivalently, at 
z = oo. 

Using these formulas and orthogonality conditions we find 

qi{t) = f\t)j^^z'-WMz,t)) (88) 
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piit) = m ^ —z'-wM^.t)) 



and 



f{t) 



dzVMz^t)) 



(89) 



(90) 



2711 z^ 

Multiplying ([89|), ( pBD by y^ and summing over all I = 0, 1, . . . oo we arrive at the final 
form of the integral equations 



, , 1 / dz V((p(z,t)) 



p{y,t) 



1 



dz vMz,t)) 



(91) 



(92) 



y J 2mz 1 - zyf{t) 

where the integration contour encircles singularities of potentials, but excludes the pole at 
Computing residues in this pole and at infinity, we arrive at algebraic form of 



yf{t) ■ 
equations 



fit) 

(i{y,t) 
p{y,t) 



vMz,t)) 



-t + Res 



;"'''H?M-'r^'^ 



^'^'KSM-'. 



Res 



z{i - zyfit)) 
vMz,t)) 



A^ - ^yfit)) 

As for the free energy, in the WKB approximation we find 

I91n/(0,x) r^^ r dz 1 



N d(t)a 



2mzp{z,t) - (j)a 



(93) 
(94) 
(95) 

(96) 



Let us now compute the second derivatives of our integral with respect to eigenvalues 
bajXb- The simplest one is the mixed derivative 



d'lnl{4>,x) 
d(padxb 



E 



1 



(Xi - (j)a){yj - Xb) , 



-E 



1 



1 



N ij 



^i 



N 



, yj - Xb , 



(97) 



N 



where averaging < A >i^ corresponds to initial measure, depending of all N x,y variables, 
rather than above averaging < A > over just one pair. 

In virtue of symmetry of the measure all terms with i ^ j contribute the same as the 
term with i = l,j = 2, and all terms with i = j contribute the same as i = j = 1, therefore 



E 



1 



1 



) = N(N -1) { ) +N' 

ij \{Xi-<pa){yj-Xb)l N \{Xi - (j)a){y2 - Xb) I j^ \{Xl - <j)a){yi - Xb) I N 



E 



1 



Xi 



N 



.y^-Xb, 



N^ 



N 



Xi - 



N 



' 1 ' 

,yi-Xb, 



(98) 
(99) 



N 



Repeating the same steps as before, when we derived (^), we find here after some 
cancellations 

d<PadXb ^l\{Xl - <Pa){yi - Xb) I n.tiA Xi - <pa /\ 2/2 " Xfe / 

where angular brackets already correspond to one dimensional averaging. 
Now we further reduce terms by applying the completeness relation 

^ / Pn{Xl)QUyi) \ / Pm{x2)Qn{y2) \ _ I Pn{x^)Qn{yi) \ ^^^^^ 

^l\ Xi-Cpa /\ y2-Xb I \iXi-(f)a)iyi-Xb)/ 

and we find after changing summation variables n = N — k,m = N + l + l 

d'lnli<p,x) _yy/ PN-kix,)QN+l+^iy^) \ / PN+i+iix2)QN^kiy2) \ ^^^3) 

d(j)adXb fro io \ Xi-(t)a / \ y2-Xb / 

So far we did not use any approximations. Now we use the WKB approximation at large N 

' PN-kixi)QN+i+iiyi)\ ^ j; dzi Zi+' 



(104) 



xi-^a / J 2inp{zi,l)-(l)a 

^N+l+l{x2)QN-k{y2)\ 



y2-Xb I 

hN+l+1 PN+l+l{x2)QN-k{y2) 



hN-k \ y2-Xb 



2mq{z2,l) -Xb 

dnnl{(j),x) 
d(f)adxb 

dzi f dz2 fiX) 



(105) 



2m J 2m (1 - f (I) z^Z2Y {p{zu 1) - 0a)(g(^2, 1) - Xb) 

The computation of the double derivative g° q^^' for a 7^ 6 goes along the same line 
untill the following point 

aMnJ(0,x) a^b , . 

o , o , = (106) 

d<pad<pb 

^ ^ / PN-kiXi)QN+l+liyi)\ I PN+l+l{x2)QN-k{y2)\ 
k=0 1=0 \ 



Xl-(j)a / \ X2- (pb 
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When we apply the WKB approximation, we have to expand P polynomials in both averages, 
unlike before, when we expanded P in first average, and Q in the second one. The result is 
therefore, slightly different 

d^\nl{(j),x) a^b j; dzi j; dz2 1 ,^^^. 



2m I 2m {zi - ^2)2(^(^1, 1) - 0a)(p(^2, 1) - ^b) 
At a = b there is an extra term, calculable by the same method, adding this term we find 
dHnI{<P,x) 



108 



N5ab I dt 



/ 2mz{q{z,t) - (paY 
dzi f dz2 1 



2m I 2m {zi - Z2y{p{zi, 1) - (f)a){p{z2, 1) - (j)b) 

However, as we found, the Sab terms in the effective Action all cancel among themselves. This 
Sab term reduces to V"{(j)) which adds to similar term from the Vandermonde determinant 
and the mass term and cancel in virtue of the classical equation. 

Let us discuss the choice of the z integration contours. Originally they enclosed the 
singularities of q{z, t) and all the zeros of denominators. After summing the geometric series 
the new double poles appear. The pole at f{l)ziZ2 = 1 must be taken outside the integration 
contours, as it arose after continuation of the geometric series from domain of convergence 
Z1Z2 — > 0. The situation with the pole at zi = Z2 is slightly more delicate. Never mind how 
small a circle you take for the ,21,-22 integrations this pole would always be located directly 
on this circle. However, the residue at this pole, say, for Z2 

^ I ^ I (109) 

2m {p{zi, 1) - (/)a) dzi {p{zi, 1) - (pb) 

would be antisymmetric with respect to a, b, as the closed loop integral of the total derivative 
vanishes. Hence, this pole should be disregarded, to preserve the symmetry of the matrix of 
the second derivatives. In other words, this is the principal value 

' ' ' 4^-^^ (no) 



{zi-Z2y 2{zi-Z2 + t0y 2{zi-Z2-t0y 
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